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Nonadiabatic pure spin pumping in zigzag graphene nanoribbons with proximity 
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By combining Floquet theory with Green’s function formalism, we present non-adiabatic quantum 
spin and charge pumping through a zigzag ferromagnetic graphene nanoribbon including a double¬ 
barriers structure driven weakly by two local ac gate voltages operating with a phase-lag. Over a 
wide range of Fermi energies, interesting quantum pumping such as i) pure spin pumping with zero 
net charge pumping, ii) pure charge pumping and iii) fully spin polarized pumping can be achieved 
by tuning and manipulating driving frequency in the non-adiabatic regime. Spin polarized pumping 
which is measurable using the current technology depends on the competition between the energy 
level spacing and driving frequency. 
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I. INTRODUCTION 

Time-periodic gate voltages applied on a conductor 
can generate a dc pumped current even at zero bias 
voltagei”— . Quantum charge pumping can be produced 
when the coherence length of electrons in mesoscopic con¬ 
ductors (e.g. graphene and nanotubes) is larger than the 
device length. On the other hand, additional to several 
micron mean-free path of electrons, graphene has also 
very weak spin-orbit coupling^ and, its spin flip length is 
so long about Ijim at room temperature^. So graphene- 
based devices have a good opportunity for spintronic 
applications^. However, weak spin-orbit in graphene re¬ 
sults in a weak electrical control over spin polarization. 
With such large spin and phase coherence length of elec¬ 
trons in graphene, it is expected that time-periodic gate 
voltages can generate significant spin pumped currents 
passing through ferromagnetic graphene^. Recently, by 
improving experimental techniques, spin splitting of elec¬ 
trons in graphene is induced by strong proximity of 
a ferromagnetic insulator layer which is deposited on 
graphene sheet^ such that its large mobility of carriers is 
preservediS. Here we apply a time-periodic electrostatic 
potential to control spin and charge transport. 

The peculiar properties of bulk graphene such as 
chiralityii induce surprisingly large pumped current 
originating from evanescent modes at the Dirac point 
which is larger than the pumped current measured in 
semiconductor s . On the other hand, controlling of 
Van-Hove singularities, ribbon size, the edge structure 
and also defects make graphene nanoribbons (GNR) and 
also nanotubes as the appropriate candidates for en¬ 
hanced pumped currenfpi^-— . Furthermore, spin current 
in monolayer graphene has been also proposed to gener¬ 
ate by application of two oscillating gate voltages through 
the adiabatic quantum pumping in some special Fermi 
energies^. However, the main point of the present work 
is the control of spin current by tuning of the driving 
frequency in non-adiabatic regime. Moreover, in this 
regime, it is possible to tune pumped current by using 
the other electrical parameters^”—. In this regime, the 
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FIG. 1: Schematic view of zigzag graphene nanoribbon 

ZGNR(M,N) with M and N atomic layers in length and 
width, respectively. The reservoirs are considered to be the 
same as the sample located at the central region. Nanorib¬ 
bons are driven by two periodic-time dependent gate voltages 
with amplitude Vac, frequency Q, and phase difference (p. Two 
large barriers (Vg) are applied in the first and last unit cells. 


possibility of inducing Floquet topological properties by 
laser illumination has recently opened a new interesting 
field^S. 

In this paper, by using Floquet Landauer 
formalisn>^i2i, we investigate a crossover from adia¬ 
batic to non-adiabatic charge and spin pumping in 
ferromagnetic graphene nanoribbons induced by two 
weakly local time-periodic gate voltages operating with a 
phase difference. Depending on the system parameters, 
it is shown that a pure spin and charge pumping emerges 
at some special driving frequencies. At low frequency 
regime, spin splitting is not remarkable such that there 
is no pure spin pumping at all, while at high frequency 
regime and also in the presence of the resonant states 
arising from the barriers, the situation for emerging pure 
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spin pumping is achievable. At the end, we present a 
discussion about practical implementations of a quantum 
pump based on graphene nanoribbons. 

This paper is organized as the following: In section 
II, we present the methodology for calculating pumped 
current by using Floquet approach in combination with 
Green’s function formalism (Floquet Landauer formal¬ 
ism). Then in section III, we present our results on charge 
and spin pumping for zigzag graphene nanoribbons. Fi¬ 
nally, the last section includes the conclusions. 

II. FORMALISM 

As shown in FigO we consider zigzag graphene 
nanoribbon in the proximity of ferromagnetic insulator 
such as EuC& and yttrium iron garnet (YIG)^ which is 
connected to two infinite electrodes. For simplicity, the 
electrodes are considered to be the same nanoribbon as 
the central system. One tt orbital is considered per each 
site for graphene system. So the single electron Hamilto¬ 
nian of the central region is described as 

Hcit) = (t) - CFh)c\^Ci„ + 

where c\^ and Cia- are the electron creation and an¬ 
nihilation operators and to is the hopping energy be¬ 
tween nearest neighbour atoms, respectively. All ener¬ 
gies are scaled in unit to- cr = ±1 refers to up and 
down-spin. Here h is the induced exchange field aris¬ 
ing from strong proximity of a ferromagnetic insulator. 
This induced exchange field can be tuned by means of 
an in-plane external electric field^S. The exchange split¬ 
ting inducing by a ferromagnetic insulator such as Euo 
in graphene has been estimated to be in the range of 
5 — 65meH— The coupling Hamiltonian which is in¬ 
dependent of time is defined as the following: Hcoupling = 

+ tRcCiC^j) + H.c., where Ilc and tnc are 
the coupling energy between the electrodes and central 
ribbon. To have quantum pumping, one should break 
left-right symmetry by application of a dynamic or static 
asymmetric spatial factor. Here, we will apply a dy¬ 
namical break of the left-right symmetry appearing as 
a phase lag between two oscillating gate voltages. There¬ 
fore, in the case of static symmetry such as 
one can still observe dc pumped current. Harmonically 
time-dependent gate voltages are applied just on the last 
and first atomic layers of the central region as shown in 

Figm 

ef (t) = eVf -h eV-'^cos{Vtt -h (fi) (2) 

where and Vi are the strength and frequency of 
the oscillating gate voltages, respectively. Moreover, let 


us apply zero gate voltage on all sites =0, except the 
first and last unit cells that include two barriers as strong 
as 1^^ = 1. Constant gate voltage causes to shift Van 
Hove singularities into the band center giving rise an en¬ 
hanced pump currenfpi^. To break the left-right symme¬ 
try, as shown in FigHJ two ac gate voltages with a phase 
lag (f = (fR — ipL are applied on the first {cpL =0) and 
last {ipR = ip) unit cells along the transport direction. We 
consider the phase difference between two ac gate volt¬ 
ages to be (/? = 7r/2 in which the maximum dc current has 
been reported in quasi-one dimensional system^i22i2^. 
Two semi-infinite electrodes are not driven by any time- 
dependent gate voltages. Because of zero source-drain 
voltage, we set on-site energies of the electrodes to be 
zero. 

To obtain dc current through graphene nanoribbon, 
we use a Floquet theory in combination with non¬ 
equilibrium Green’s function formalism. It should be 
noted that the electron-electron interaction is not consid¬ 
ered in this calculation. At zero temperature, the aver¬ 
aged pump current over one period of the ac gate voltage 
for spin a is written as the following^ i^^i^^ . 

C = iEnI-:.deT^i^) (3) 

where 

T:{e) = TRcAe)-TEnAe) 

and Er IS the Fermi energy. The transmission probability 
of carriers coming from the left to the right electrode 
T'^l a(^) assisted by the absorption (n > 0) or emission 
(n < 0) of I n I photons with an incident energy of carriers 
5 and spin a can be written in terms of Floquet Green’s 
function G^{s) as the following. 

TRlA^) = Tr[rRje)GRLA^)rioi^)GiLAs)] ( 4 ) 

where Floquet Green’s fucntion is calculated in based 
on the Floquet states as ^{e) = {R^n\{el — 

and Floquet Hamiltonian is defined as = 
Hc{t) — The Floquet Green’s function is a Fourier 

transformation of retarded Green’s function as the fol¬ 
lowing: 

oo 

GRL,At,£)= E 

n=—OO 

Taking Fourier transform of the equation of motion 
governing on the retarded Green’s function results in an 
equation of motion followed by the Floquet Green’s func¬ 
tion. The escape rates of electrons to the electrodes 
through elastic or inelastic channels are defined as the 
following. 

ra,n(e) = + nHuj) - + nhuj)] 



3 





FIG. 2: Pumped current in terms of driving frequency at the 
band center = 0 for ZGNR(12,4)- Blue solid line (red dash- 
dot line) shows charge (spin) pumped current. In the lower 
panels, some of frequency’s ranges are magnified to present 
some special frequencies. 


where the index n shows the number of photons as¬ 
sisted for inelastic transport and Hq, {e) is the Floqnet self 
energy arising from the electrodes a = L^R. In the cal¬ 
culations, the convergence of the transmission functions 
determines the number of required Floqnet channels n. 
In the case of non-interacting systems, the above for¬ 
malism is equivalent to non-equilibrium Green’s function 
formalism (Keldysh formalism)^^. 


III. RESULTS 

A. Driving Frequency Dependence of Pumped 
Current 

We present in FigE] the results for spin and charge 
pumped current through zigzag graphene nanoribbons 
as a function of driving frequency Q at the band cen¬ 
ter. The charge and spin pumped current are defined as 
Jc = /t + h respectively. and 

refer to the up and down spin currents. For very low 
frequencies (adiabatic regime), there is a linear behav¬ 
ior of charge pumped current in terms of Q which is in 
compatible with the results reported beforei^. This lin¬ 
ear behavior which is shown in Fig.(|2lb), is a character 
of the adiabatic pumpingi. However, close to the Van 
Hov singularities, the scale of pumped current with the 
frequency is improvedi^. 

The characteristic energy scale in this system is de¬ 
fined as Umax * ^ in which Umax IS the maximum number 
of sidebands excited by the driving gate voltages. The 
adiabatic regime occurs when this energy scale is small 
enough in compared to the level spacing of the resonant 
states. On the other hand, Umax depends on the ratio 
of Vac/^- It means that the adiabatic regime occurs at 


lower frequencies when the strength of pumping is large. 
Furthermore, longer or wider ribbons result in smaller 
level spacing. Therefore, larger pumps have smaller fre¬ 
quency ceiling for working in the adiabatic regime. The 
adiabatic regime is indicated by the point 1 in Fig.(j2lb) 
which corresponds to the pumping frequency as small as 
= O.Olto* Although the level spacing is energy depen¬ 
dent, in our system, one can estimate it to be about 0.2 
(to) which is larger than U = 0.01 (to). In the adiabatic 
regime, spin pumping is small, however there is no pure 
spin pumping in this regime. 

At high frequency regime, spin and charge pumped 
current show some extremum points at the special val¬ 
ues of the driving frequencies which are originated from 
the resonance of the Floqnet states with those atomic 
levels captured between two barriers. Therefore one can 
control the direction and amplitude of the spin current 
by tuning the driving frequency. At frequencies that are 
multiples of the level spacing, by mediated of photon- 
assisted transitions, atomic levels are mixed with each 
other. Therefore, an enhanced pumped current emerges 
in this situation. However, level spacing is not uniform 
along the whole spectrum, so an interplay between driv¬ 
ing frequency and level spacing is more complicated. 

In the non-adiabatic regime, spin-polarized pumped 
current is observable along the whole range of driving fre¬ 
quencies shown in Fig.(|2]). However, full spin polarized 
pumping are categorized in the following three interesting 
cases: i) in the case of Ic = Is^ there exists pure up-spin 
polarized current. So we have 7 ^ 0, /^ = 0. This case is 
marked by the point 2 in Fig.(| 2 ]) referring to U 2 . If charge 
and spin pumped current are positive/negative, up-spin 
current is pumped from the left/right to the right/left 
electrode, ii) The case of Jc = 0 which emerges at the 
frequency U 3 , results in a fully spin pumping with zero 
net charge pumping. The pumped currents for up and 
down spins are in the opposite directions 

In fact, quantum interference is responsible for such 
interesting cases. In the adiabatic regime, pumped cur¬ 
rent is an odd function of the phase difference between 
oscillating gate voltages. However, in the non-adiabatic 
regime, an additional contribution to the phase difference 
is originating from the spatial phase difference of the Flo- 
quet sidebands2£. This phase difference which is absent 
at low frequency regime, is responsible for the sign re¬ 
versal of the pumped current at consecutive peaks shown 
in Fig El Photon-assisted Quantum interference through 
different Floqnet side-bands can result in a special case; 
two spin-polarized pumped current which flow in the op¬ 
posite directions. So in this case, there is no net charge 
pumping while it is purely spin-polarized pumping in the 
opposite directions. 

hi) Correspondingly, the case of Ic = —Is also results 
in a fully down-spin polarized currents I^ = 0 ,/^^ 7 ^ 0 . 
If/c > 0, then we have I^ > 0. This frequency is in¬ 
dicated by U 4 . In other cases, the pumped current is 
partially spin polarized. Spin polarization emerges when 
the density of state for up and down spins shift against 
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FIG. 3: Pumped current in terms of Fermi energy 

for ZGNR(12,4) at the specified driving frequencies of 
Qi, ^ 2 , which have been marked in Fig[2 
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FIG. 4: Net pumped transmission in terms of Fermi en¬ 
ergy for ZGNR(12,4) at the specified driving frequencies of 
, Q 2 , which have been marked in FiglJ] The norm of the 
left-right green’s function for single photon-assisted transmis¬ 
sion is plotted in part (b) for the low frequency regime Qi. 


each other, so at a given frequency, level spacing for up 
spins differs from one for down spins. As a consequence, 
pumped current for up and down spin would be differ¬ 
ent in its value and also pumping direction. Moreover 
the effect of higher orders of the orbital overlaps on the 
pumped current is investigated in the Appendix. 


B. Fermi Energy Dependence of Pumped Current 

Let us concentrate on the Fermi energy dependence of 
the pumped current at these four special driving frequen¬ 
cies shown in Fig|3l Panel (a) in Fig l3] shows spin-up and 
down pumped current as a function of Fermi energy for 
small pumping frequency = 0.01 to- This frequency is 
smaller than two consecutive energy levels of the cen¬ 
tral region. The pumped current contains successive 
small peaks which are modulated on a pumped currenij^^. 
These peaks in the pumped current are originated from 
the peak-antipeak structure^ of the net pumped trans¬ 
mission T^) shown in Fig.dHa). In this regime, there 
is no mixing between system levels, therefore spin po¬ 
larization is small. However, there exists unidirectional 
charge pumping along the whole range of Fermi energies. 

Panels (b), (c) and (d) in Fig IS indicate fully spin po¬ 
larized pumped current in terms of Fermi energy for the 
frequencies 112,^3 5 ^ 4 - Apparently, because of two bar¬ 
riers applied on the first and last unit cells, electron-hole 
symmetry breaks. Furthermore, the pumped current in¬ 
creases around E" = — 1 which is in resonance with the 
barriers. In panel (b) and at (^ 2 , down-spin current is 
blocked and up-spin is fully pumped from the right to the 
left electrode. At and in the panel (c), we present an 
interesting case of pure spin current with zero net charge 


pumping over a wide range of Fermi energy around the 
band center. Finally, at the frequency ( 24 , down-spin is 
fully pumped, while up-spin is still blocked. There is 
some structure of jumps and plateaus in panels (b) to 
(d) which follow the resonant states shown in the net 
pumped transmission curve (Figl4l3,c). The behavior of 
pumped current in ZGNRs is similar to the pumped cur¬ 
rent passing through one-dimensional chain^. 

Fig.© shows the net pumped transmission curve 
(which is defined as as a function of energy for 

different driving frequencies. For low driving frequency 
(2i, peaks and anti-peaks in the net transmission curve 
emerge around the resonant states located between two 
barriers (Fig.(j4la)). At both energies of peak and anti¬ 
peak, we observe large quantum pumping and hence con¬ 
structive interference. On each peak-antipeak pattern, at 
lower energies, the left to the right pumping is dominant 
while at higher energies, pumping direction will be from 
the right to the left side. The following adiabatic picture 
may be useful to understand the reason for such behav¬ 
ior: At half period of oscillating potential, the left barrier 
has a lower gate voltage in compared to the right barrier. 
So electrons flow from the right reservoir into the well. 
In the other half of period, the right barrier gets lower 
which leads to a current flow from the well into the right 
reservoir. If we turn off the oscillating gates, transmission 
would have peaks at the resonant states as the original 
transporting channels. 

At low frequency regime, there is no mixing between 
system levels, so transport occurs mainly through reso¬ 
nant and also one Floquet side state. However, at high 
frequencies, the system levels are mixed with each other 
arising from the driven potential. As a result, large peaks 
are appeared in the net pumped transmission which are 
observed in Figs.(|llb,c). 
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C. Practical Implementations 

In reality and in the absence of coulomb blockade, stray 
capacitances between the gates and electron reservoirs 
may affect adiabatic quantum pumping in the work done 
by Switkes^. However, some features of the observed 
dc voltage are in agreement with the theoretical calcu¬ 
lations based on the scattering matrix—. To have con¬ 
clusive spin pumping, one of the powerful techniques is 
precessing magnetization at ambient temperature^. Fur¬ 
thermore, the other technique is the use of the ac Joseph- 
son effect in a hybrid normal-superconducting systent^or 
charge pumping in an unbiased InAs nanowire embedded 
on a superconductor—. For graphene nanoribbons with 
smooth edges, in the adiabatic regime, it is expected to 
observe some general feature of quantum pumps at low 
temperatures. Furthermore, we have shown before^ that 
in graphene nanoribbons, under application of a weak 
source-drain voltage, screening effects are weak. There¬ 
fore, in a weak driving gate voltages, accumulation and 
depletion of charge would be small. On the other hand, 
homogeneity of the junctions reduce stray capacitances 
between oscillating gates and electron reservoirs. Here 
we assume that the electrodes are similar to the pump¬ 
ing region for having transparent contacts in compared 
to the normal pumps. 

In reality, the width of nanoribbons are much wider 
than what we have considered in this paper. Wider and 
longer nanoribbons have smaller energy level spacing. So 
it is expected that for wider and longer nanoribbons, 
more sign reversal happens for the pumped current curve 
in terms of the driving frequency. For graphene sheet, 
there exists pure spin pumping at special energies in the 
adiabatic regime^. Pure spin pumping is also expected 
to observe in wider graphene nanoribbons. 

During fabrication processes, the ribbons are full of 
impurity or disorder. So we propose to investigate the 
effect of disorder and impurity locations on the pure spin 
and charge pumping through graphene nanostructures. 


IV. CONCLUSION 

In conclusion, we study charge and spin pumping in 
ferromagnetic graphene nanoribbons including a double¬ 
barrier structure driven weakly by two local ac gate volt¬ 
ages with a phase-lag. With new improvements to induce 
exchange magnetic field in grapheneiS, the main point of 
this paper is that charge and spin pumping can be well 
controlled by tuning the driving frequency in the non- 
adiabatic regime. Furthermore, fully spin pumping with 
favorite direction is achievable at the special frequency of 
ac gate voltages. It seems that at low frequency regime, 
single photon-assisted tunneling is dominant, while at 
higher frequencies, mixing of energy levels arising from 
the driving field leads to constructive interference. 



Driving Frequency 


FIG. 5: Spin pumped current (Is) as a function of driving 
frequency for ZGNR(20,4) at the band center for two cases: a) 
by considering of the nearest neighbor overlaps (NN-GNR), 
b) the second and third nearest neighbor hopping energies 
(3NN-GNR). 
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VI. APPENDIX 
A. Higher order orbital overlaps 

In the tight binding model shown in equation [H hamil- 
tonian contains nearest neighbor interaction. However, 
ab initio band structures can be extracted by consider¬ 
ing of the second and third nearest neighbor overlaps 
in Hamiltonian. To investigate the effect of higher or¬ 
der overlaps, the pumped current is calculated by using 
tight-binding parameters which are proposed by Ref. (Isbl ) . 
The distances between the first, second and third nearest 
neighbors are given by ai = 1.42A,a 2 = \/3ai,a3 = 2ai. 
The values of nearest-neighbor hopping parameters are 
calculated by comparing and fitting the first principle 
and tight-binding band structures for different types of 
graphene^ (2D graphene, graphene nanoribbons and 
nanotubes ). The second and third hopping parame¬ 
ters are considered as 0.12eV and 0.068eV if the nearest- 
neighbor hopping to is considered to be as 2.78eV— . 

Fig l5] shows spin pumped current (Is) in terms of driv¬ 
ing frequency for the two cases: Considering of a)the 
nearest neighbor overlap (NN-GNR), b) the second and 












third nearest neighbor overlap (3NN-GNR). As it is the 3NN-GNR case in compared to the NN-GNR case, 
shown, spin pumped current shows more extremum for 
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